Field— and pressure— induced magnetic quantum phase transitions in TlCuCls 
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Thallium copper chloride is a quantum spin liquid of S = 1/2 Cu'^"'' dimers. Interdimer superex- 
change interactions give a three-dimensional magnon dispersion and a spin gap significantly smaller 
than the dimer coupling. This gap is closed by an applied hydrostatic pressure of approximately 
2kbar or by a magnetic field of 5.6T, offering a unique opportunity to explore the both types of 
quantum phase transition and their associated critical phenomena. We use a bond-operator for- 
mulation to obtain a continuous description of all disordered and ordered phases, and thus of the 
transitions separating these. Both pressure- and field-induced transitions may be considered as the 
Bose-Einstein condensation of triplet magnon excitations, and the respective phases of staggered 
magnetic order as linear combinations of dimer singlet and triplet modes. We focus on the evolution 
with applied pressure and field of the magnetic excitations in each phase, and in particular on the 
gapless (Goldstone) modes in the ordered regimes which correspond to phase fluctuations of the 
ordered moment. The bond-operator description yields a good account of the magnetization curves 
and of magnon dispersion relations observed by inelastic neutron scattering under applied fields, 
and a variety of experimental predictions for pressure-dependent measurements. 

PACS numbers: 75.10.Jm, 75.40.Cx, 75.40.Gb 



I. INTRODUCTION 

Thallium copper chloride^'^'^ presents an insulating, 
quantum magnetic system of dimerized 5 = 1/2 Cu^+ 
ions. Inelastic neutron scattering (INS) measurements 
of the elementary magnon excitations'*'-^ reveal a strong 
dispersion in all three spatial dimensions indicative of 
significant interdimer interactions. The dispersion min- 
imum gives a spin gap Ag = O.TmeV, which is signif- 
icantly smaller than the antiferromagnetic (AF) dimer 
superexchange parameter J « 5meV. The correspond- 
ing critical field, He — 5.6T, makes TlCuCls one of the 
few known inorganic systems in which the gap may be 
closed by application of laboratory magnetic fieldsi^ Neu- 
tron diffraction measurements at fields H > He revealed 
that a field-induced AF order in the plane normal to 
the applied field appears simultaneously with the uni- 
form moment 1^ Recent INS measurements of the magnon 
spectra in finite fields,^ including those exceeding Hc,^ 
have provided dynamical information concerning the el- 
ementary excitations, in particular the linear Goldstone 
moderi in the phase of field-induced magnetic order. 

TlCuCla (Fig. 1) is one member of a group of related 
compounds. The potassium analog KCuCl j*i^i''i*°i-^^i^^i^^ 
is similarly dimerized, but has significantly weaker inter- 
dimer couplings, resulting in a large spin gap of 2.6meV. 
A further material in the same class, NII4CUCI3, has no 
spin gap and exhibits magnetic order with a very small 
moment, but also shows a complicated low-temperature 
structure which gives rise to magnetization plateaus only 
at 1/4 and 3/4 of the saturation value. While the ap- 
parent increase of interdimer couplings with anion size 
may suggest a contribution of the anion to superexchange 
processes, it should be noted that the physical origin of 



the properties of NII4CUCI3 may be rather different from 
the other members.'''^ Turning from chemical to physical 
pressure, Tanaka et al}^ found by magnetization mea- 
surements under hydrostatic pressure that TlCuCla has 
a pressure-induced magnetically ordered phase, with a 
very small critical pressure for the onset of magnetic or- 
der, Pc = 2kbar. Oosawa et ali^ have shown very re- 
cently by elastic neutron scattering measurements under 
a pressure of 1.48GPa that the pressure-induced ordered 
phase has a strong staggered moment (60% of the satu- 
ration value), again reflecting the low value of Pc- The 
magnetic Bragg reflections are found at reciprocal lattice 
points Q = (0,0, 27r) (following the notation of Ref. 0), 
as in the field-induced ordered phase of TlCuCla. The 
aim of the present work is to compare and contrast the 
field- and pressure-induced ordered phases of the sys- 
tem, and to provide a complete description of the static 
magnetization and dynamical excitations at all fields and 
pressures. 

The issue of coupled spin liquids and closing of the 
spin gap due to higher-dimensional superexchange inter- 
actions has been investigated in a wide variety of sys- 
tems, many of which were found to be rather close to the 
quantum phase transition (QPT) to an ordered state. 
The Haldane chains CsNiCla and RbNiClg:*^. presented 
an early example where both static and dynamic prop- 
erties could be measured. The coupled, two-leg ladder 
system LaCu02.5*^ was found to have sufficiently strong 
interladder coupling that magnetic order could set in, 
whereas the coupled plaquette compound CaV409^— re- 
tains a robust spin gap. While numerical simulations of 
appropriate models have compared the thermodynamic 
properties of the phases on both sides of the transition^Si 
theoretical studies of these systems have concentrated 
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in particular^^i^^i^'*i^^i^^i^'''i'^^ on the identification of a 
low-lying but massive longitudinal mode on the ordered 
side of the transition which corresponds to an ampli- 
tude fluctuation of the ordered moment. Experimental 
work continues to refine the observations of this type of 
excitationi^S 

The first system to be discovered in which the spin 
gap could be closed by application of a magnetic field 
was the organometallic compound (C5Hi2N2)2Cu2Cl4i^ 
This material, originally believed to have a two-leg lad- 
der geometry, undergoes a quantum phase transition 
(QPT) at He — 6.6T, and has been investigated quite 
extensiveljiSL^SiS to characterize the thermodynamic 
properties in the quantum critical regime. The magnetic 
interactions within this system are now thought to have 
a significantly more complex structure?^ whose 3d na- 
ture mandates a reinterpretation of earlier observations. 
A further recent discovery, (C5Hi2N2)2CuBr4(2^ is cur- 
rently believed to be a quasi-ld ladder system, and also 
shows an almost identical lower critical field. He = 6.6T. 
A considerable body of theoretical ■work2Si2L^2i224£ has 
analyzed the properties of this system by a variety of 
techniques. A particular focus of recent studies in this 
direction, the majority motivated by experiments on 
TlCuCls, has been the possibility of describing the or- 
dered phase as a field- or pressure-driven Bose-Einstein 
condensation (EEC) of the triplet magnon excitations of 
the spin liquidj^ii^SiiSiii which become massless at the 
QPT. In this formulation the critical exponents of ther- 
modynamic properties such as the magnetization in the 
ordered phase, are specified as a function of dimensional- 
ity and magnon dispersion, and allow in addition a char- 
acterization of the quantum critical regime at finite tem- 
perature. However, questions remain concerning of the 
appropriateness and universality of this description, and 
we will address these below. 

Here we work exclusively at zero temperature to con- 
sider the field- and pressure-induced QPTs in TlCuCla 
using the bond-operator technique. This approach is 
naturally suited to the spin-liquid regime of dimerized 
systems^ and has previously been developed both for 
coupled ladders"'^ and for describing the neighboring 
magnetically ordered phases2Li2ilS on an equivalent foot- 
ing. This analytical framework is thus uniquely applica- 
ble for a consistent discussion of static properties, and 
in particular of dynamical excitations, on both sides of a 
QPT. In Sec. II we consider the structure of TlCuCla 
and deduce the parameters of the effective model for 
the elementary magnon excitations in zero field. Sec. Ill 
contains the bond-operator description of the quantum 
phases and phase transitions in the TlCuCls system as 
a function of the field, and as such constitutes the ex- 
tended version of our previous results for magnetizations 
and magnon excitationsi^ In Sec. IV we apply the bond- 
operator formulation to the pressure-induced QPT, ex- 
tending the treatment of Refs. and to discuss in 
detail the phase and amplitude modes in TlCuCla, the 
effects of simultaneous pressure and field, and finally the 




FIG. 1; Structure of TlCuCls: small circles represent Cl~ 
ions, medium-sized circles Cu^^ ions, and large circles Tl"*" 
ions. 



magnetization at finite pressure. Sec. V presents a sum- 
mary and conclusions. 

II. MAGNETIC BAND STRUCTURE 

TlCuCla and KCuCls appear to have the structure of 
a quasi-ld system of coupled spin chains. However, INS 
experiments on both systems at zero field reveal an el- 
ementary triplet excitation with a spin gap and strong 
dispersion in all three reciprocal-space directions. In this 
section we using a model of coupled dimers to provide a 
qualitative justification for the observed results in terms 
of superexchange interactions within the crystal struc- 
ture. Although the dimer model is different from those 
of Refs. y, |5j, and ll^ the results are perforce in close 
agreement. 

A. Exchange Pathways 

The magnetic couplings between the spins, which are 
due to unpaired electrons in the Cu d,j.2_y2 orbitals, are 
mediated by superexchange through the p orbitals of the 
Cl~ ions. Quantitative calculations of superexchange in- 
teractions, particularly for extended pathways, remain 
unfortunately beyond the scope of current understand- 
ing and computer power. Here we employ qualitative 
considerations of geometry and orbital overlap to provide 
plausibility arguments for the relevant model parameters. 

Each Cu ion can be considered to be octahedrally co- 
ordinated by CI (Fig. 1). The Cu dx2_y2 orbitals lie in 
the basal planes of the octahedra, and the most impor- 
tant superexchange interaction, J, is that between ions 
in two edge-sharing octahedra. Such a pathway requires 
exchange through quasi-orthogonal orbitals on CI, and if 
the environment is truly cubic Hund's Rule dictates that 
the coupling is weakly ferromagnetic (FM). However, for 
a structure distorted away from this high symmetry, as in 
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FIG. 2: Schematic representation of relevant interdimer cou- 
plings in XCuCla; (a) a-c plane; (b) b-c plane. 



TlCuCls, AF exchange terms are expected to be stronger, 
although the edge-sharing configuration remains a weak 
superexchange path on the scale of the linear bonds in 
cuprates. Further superexchange terms in the chains (left 
side of Fig. 1) will of necessity involve pathways coupling 
to the Cu (i3^2_^2 orbitals, and thus will be considerably 
weaker than J, while interchain pathways, which have the 
form Cu-Cl-(T1)-Cl-Cu (Fig. 1), should also be small. 
From this consideration one expects a system composed 
of dimers, provided by the edge-sharing octahedra, with 
weak interdimer coupling in all directions. 

Nearest-neighbor couplings in the chain-like stack of 
dimer units (oriented along a) are also of the edge- 
sharing type between octahedra, where the edge is from 
basal plane to apex. Again the possibility arises from the 
geometry that these weak bonds may be FM. Viewing the 
dimers as ladder rungs, the effective ladder leg coupling 
is denoted by Ji , and the asymmetrical cross-coupling by 
J'l [Fig. 2(a)]. Finally, the structure of the ladder legs is 
such that a next-neighbor Cu-Cl-Cl-Cu coupling, sim- 
ilar to that proposed in CuGeOs, may not be excluded, 
and this is denoted as J'{ . A significant interchain cou- 
pling is expected only where the Cu-Cl-Cl-Cu pathway 
is rather direct, with bond angles not deviating widely 
from 180*^. The dimers are oriented primarily along the 
c axis (Fig. 1), where by inspection any overlap is not es- 
pecially direct and may involve parallel paths. Because 
the dimers are also tilted significantly along a, it is pos- 
sible that dimers with a separation of one or two unit 
cells in this direction are coupled; let J2 denote a dimer 
coupling along c only, a coupling along a + c and J2 
a coupling along 2a + c [Fig. 2(a)]. Finally, the shortest 
interchain distance is that between the two inequivalent 
chains in the unit cell, namely \c (Fig. 1). Here the 
directionality is evident from the fact that the pathway 
goes through the CI ions pointing maximally up or down 
(relative to the a-axis) from the dimers, indicating that 
this coupling, J3, should correspond to the interdimer 
separations 770-1-^6-1- ^r/c, where 77 = ±1. The possibil- 
ity remains that the couplings of the spin on one dimer 
to either of those on a neighboring dimer are finite, re- 
quiring a final superexchange parameter, J3 [Fig. 2(b)]; 



inspection of bond lengths and angles (Fig. 1) suggests 
that Jg < J3. 

B. Model 

Having argued that the properties of the system in 
zero field, including the spin gap, are determined largely 
by the dominant dimer coupling J, one may construct 
a model with this ground state using the bond-operator 
technique.'*^ In this representation the S = 1/2 spin de- 
grees of freedom on each dimer, or ladder rung, are ex- 
pressed by the bond operators 

\s) = s^\Q) = ^^{\U)-\i])), 

lU = 4l0) = -^(ITT>-|U», 

\ty) - 4lo> = 7i(ITT) + lii)), (1) 

%) = 4lo)-;^(in) + iiT», 

where the arrows in each state denote the direction of the 
two spins, and all the states created in Eq. |Q have total 
(dimer) Sz = 0. From the action of the spin operators 
S; and Sf. (denoting left and right sites of the dimer) on 
these states one may deduce"*^ the correspondence 

5"^ — iti(s^t(j -|- i]js) — ieap-ftjjt^, (2) 

where a,P,"f = x,y,z. The spin commutation relations 
for S/,r- are recovered if the singlet and triplet operators 
have bosonic statistics. However, the bond operators are 
constrained by the number of physical states available on 
each dimer, one singlet or one of three triplets, to obey 
the condition 

4s^+ E iUhc = l, (3) 
a—x,y,z 

which specifies a hard-core bosonic nature. 

With a view to the analysis to be presented in the sub- 
sequent sections of the system under a finite magnetic 
field, we introduce here the appropriate bond-operator 
description. The operators in Eq. are not the eigen- 
states of an external magnetic field, and so we transform 
these to the operators^ 

~ "71^^+ '^^-^^ 

U = to, (4) 
which create the triplet states 

tl\o) = -ITT), 

4lo) = ^(ITi) + ITi)), 

^^|o> - lii). (5) 

In this basis the magnetic-field term for a dimer has the 
diagonal representation 

^h- {S1+S2) = iKXtlty-tU,) = h{tlt+-tU^), (6) 
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in which h = g^^H , ensuring that the operators tj 
reproduce the energy levels of the field eigenstates — 
0,±1. Further, the transformation conserves particle 
number, and the local constraint retains the form 

Q = +,0,- 

FoUowing Ref. 46, transformation to the bond-operator 
representation yields 

H = H„ + Hi+H[+H'^ + H2+H^ + H^ + H3+H'3, (8) 

where 

Ho = '^[eisslsi + ^ eiatla'tta + f^i]' 

T-Ci = \ Ji \Hst{hi + Qi) + 'Htt{hi + ai)] , 

i 

K = i"f(5I[-^>^t(«,i + ai)+Htt(i,« + ai)], 

i 
i 

'H2 = \ J2 ^ [~Ust{h 1 + 0.2)+ Uttii, i + 0.2)] , 

i 

K = iJ^^[-H,t(*,i + a2)+^tt(*,« + a2)], 

i 
i 

^3 = \J?.'^\Hst{i,i + a'i±)+Utt{i,i + a3±)], 

i,± 

^3 = iJ3l][-^st(i,i + a3±)+Htt(j,« + a3±)]. 

In Tio, Cis and 6^0, denote the singlet and triplet ener- 
gies in the presence of an external magnetic field directed 
along the z-axis, and /^^ is a chemical potential emerg- 
ing from the Lagrange multiplier introduced to enforce 
the constraint Q on each dimer. The bond Hamiltonian 
Hst{i,j) is specified by 

Wst(i, j) - J2(*ltJC.4s^ + tlt%S,S, + H.C.), (9) 
a 

where a = 0, — and a = — , 0, and describes singlet- 
triplet interaction processes: the first two terms in Eq. @ 
correspond respectively to triplet propagation and triplet 
pair creation. The triplet-triplet interaction term is 
given by 

Wtt(*, j) - [tlUoitl+t,+ + + H.c] 

-[UotMA+t]- + tU]+) + H.c] (10) 

and encodes the possible triplet scattering processes. As 
in previous analyses2Li& we have chosen to neglect terms 



of the form tl^t'-^^tj^Si: these three-triplet interactions 
vanish identically for a centrosymmetric system^^*^ as 
is the case for TlCuCls in the disordered phase, and in 
the ordered phase their contributions may be assumed to 
be small. Finally, specifies the bond vectors 

di = d, 

ai = 2d, 

d2 = c, 

0,2 = d + c, 

= 2d + c, 

dsi = ±d-f iS± ic, (11) 

for the superexchange interactions of Fig. 2, where d, 6, 
and c are unit vectors for the crystal axes. The prefactor 
of each term in Eq. © is determined by the number of 
bonds between dimers, and the sign of Hst by whether 
the bond is between like {II, rr) or unlike spins [Ir). 

C. Conventional bond— operator formulation 

At zero field, as for all magnetic fields below the critical 
field He, the system is in the quantum disordered regime 
with a spin gap between the singlet and triplet states on 
each dimer. This situation, for which the bond-operator 
techique is most directly applicable, is represented by ne- 
glecting the dynamics of the singlet operator and replac- 
ing Si everywhere by a c-number, 1^, corresponding to 
a condensate of dimer singlets. Two approaches for the 
treatment of the local constraint {Tj) in the disordered 
phase exist in the literature, and we begin with a brief 
discussion of the conventional bond-operator theory,"*^ in 
which the chemical potential is retained explicitly. Pro- 
ceeding within a mean-field approximation, the opera- 
tors Si and the site-dependent chemical potentials pLi are 
replaced by uniform, global average values {si) — s and 
fii = iJ,. These parameters are then determined self- 
consistently from a minimization of the total energy at 
fixed total boson number. The dimer energies eis and eia 
in Tip take the values 

eia = jJ-ah~fi, (12) 

and the self-consistent solution follows exactly the treat- 
ment of Refs. l4^ and 27, which also present the situation 
for a system with two dimers per unit cell. 

Quantum fiuctuations about this dimer-singlet ground 
state are contained in the triplet operators in Eqs. © 
and and in reciprocal space are given to quadratic 
order by the Hamiltonian 

H = N{~ps^ - ^is^ + n) (13) 

+ X! ) X! ''katka^ka + i(^fc*fcO^-fcO + H.C.) 
+ ^kitk+t-k- + tk-t-k+ + H.C.)} , 
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where N is the number of dimers in the system. The 
dispersive part of the magnon energies for the TlCuCla 
system is contained in the expressions 

efca = jJ - ^J. + Ak- ah, (14) 
Afc =^ ^'^{(^1 ^ ^J'l) cos kx + J" cos 2kx 

- i [ J2 cos + J2 cos(fca; + kz) 

+ J2 cos(2fca; + k^)] 

+ ( J3 - J3) cos (fc^ + ifc^) cos ^ky}, (15) 

where k^ takes values in the interval — 27r < A:^ < 27r 
in order to describe the two dispersion branches corre- 
sponding to the two-sublattice system. The Brillouin 
zone lies between — tt and tt (in units of the inverse lat- 
tice spacing) in the other two directions. Diagonalization 
of the Hamiltonian H13|) by Bogoliubov transformation 
gives three modes with dispersion relations + h, i?k, 
and E'k — h, where 

Ek = /(iJW^+A^f^ (16) 

is the dispersion relation for the field-free system. In 
the disordered regime the triplet modes do not change 
the form of their dispersion, and are merely split by the 
magnetic field due to the Zeeman interaction. 

D. Holstein— Primakoff expansion 

The alternative means of treating the global constraint 
is to replace the singlet operator by using the constraint 
itself, 



y i,Q=+,o,- 

which leads to the simultaneous elimination of the vari- 
able Substitution of Eq. (|17|l into the Hamiltonian, 
followed by truncation at quadratic order after the ex- 
pansion of the square root, is equivalent to taking s = 1 
and /.t = — |J in the formulation of the previous subsec- 
tion, and at higher fields (beyond He) to specifying that 
the condensate amplitude remains unity. With these val- 
ues of s and the results of this approximation have 
exactly the form of Eqs. 

The Holstein-Primakoff approximation thus neglects 
quantum fluctuations to triplet dimers within the sin- 
glet condensate, and as such constitutes a stronger ap- 
proximation than the approach of the preceding subsec- 
tion, in which these fluctuations are included in the self- 
consistent determination of s. The values of s obtained 
for the disordered phase in TlCuCla and KCuCla are 
shown in Table I, and demonstrate that the triplet fluc- 
tuations are always small. The maximal error incurred 
is of order 3% in this regime, while at saturation, where 
quantum fluctuations are completely suppressed, the con- 
densate amplitude is in any case forced to unity. In the 





TlCuCb 


KCuCIg 


J/meV 
Ji/meV 
J 2 /meV 
J 3 1 me V 


5.50 ( 4.87 ) 
-0.43 (-0.41) 
3.16 ( 3.31 ) 
0.91 ( 1.00 ) 


A C\C\ / A '\ \ 

4.22 ( 4.16 ) 
-0.42 (-0.48) 
0.79 ( 0.75 ) 
0.70 ( 0.68 ) 


Ao/meV 


0.70 ( 0.71 ) 


2.60 ( 2.65 ) 


s 


1.00 ( 0.97 ) 
0.75 ( 0.86 ) 


1.00 ( 1.00 ) 
0.75 ( 0.77 ) 



TABLE I: Dimer superexchange parameter, J, interdimer su- 
perexchange parameters, Ji, gap, Ao, singlet condensation 
parameter, s, and chemical potential, \i, determined from 
the bond-operator analysis of TlCuCla and KCuCla . Values 
to the left in each column are obtained from the Holstein- 
Primakoff formulation of Sec. IID, while values in parenthesis 
are obtained from the conventional bond-operator theory of 
Sec. lie. Ji denotes the combination 2Ji — J'\ and J3 the 
combination Jz — J^- 

light of this observation, truncation at quadratic order, 
which is equivalent to the neglect of magnon-magnon in- 
teractions in both disordered and ordered phases, is also 
not a source of significant error. The superexchange pa- 
rameters fitted within each of the two approaches (Ta- 
bic I) allow a further comparison, and are also in rather 
close agreement, as discussed in detail in the follow- 
ing subsection. Thus the Holstein-Primakoff treatment 
yields semi-quantitative accuracy for 3d systems such as 
TlCuCls and KCuCla at low temperatures, and we will 
use this level of approximation in all of the analysis to 
follow. 



E. Superexchange parameters 

The zero-field dispersion relations H16|l deduced from 
the model in Fig. 2 may be compared directly with the 
experimental results for TlCuCla and KCuCla, which we 
take respectively from Refs. and 0. The most char- 
acteristic features of the dispersion data for both sys- 
tems are that local minima appear at both the and tt 
points in both kx- and /c^-directions, and that the rel- 
ative depths of these minima in k^ (fc^) are exchanged 
between and tt in kz (kx)- Without attempting a sys- 
tematic optimization procedure it is possible to find pa- 
rameters (Table I) giving a very good account of the 
data. Fig. OJa) shows the calculated magnon dispersion 
for TlCuClg with the data of Ref. and Fi^Efb) the 
dispersion for KCuCls with the data of Ref. Il2 . 

Under the assumption that the ±a-f -ibiic coupling is 
indeed dominated by J3, the signs of the superexchange 
interactions are consistent with the expectation that all 
intersite parameters {J} in Fig. 2 are AF. The values 
for KCuCla are consistent with the results of Refs. 0, 
and lTsL with discrepancies arising from the exact for- 
mulation of the dimer description. The leading differ- 
ence between TlCuCls and KCuCla lies in the value of 
the interladder coupling , which is larger for TlCuCls, 
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B C D A F C 

FIG. 3: (a) Magnon dispersion £k for TlCuCls. The points 
represent experimental results measured at T = 1.5K, taken 
from Ref. 0, and the solid lines the fit using the parameters 
obtained within the Holstein-Primakoff formulation listed in 
Table I. (b) Magnon dispersion for KCuCla. Experimental 
data, measured at 5K, is taken from Ref. Il2l and the theo- 
retical fit from the Ifolstein-Primakoff parameters of Table I. 
Following Ref. 0, points in reciprocal space are denoted by B 
= (0,27r,27r), C = (0,0,27r), D = (0,0,0), A = (7r,0,0), and F 
= (7r,0,27r). 



and is responsible for the contrast between wide-band, 
small-gap and narrow-band, large-gap magnon disper- 
sions exemplified by the two systems. We now discuss 
the dispersion relations and superexchange parameters 
in further detail. 

The dominant contribution to the overall dispersion 
curves conies from J3, specifically its cosfc^; cos ^fc^ form. 
Although the magnitude of this coupling is not large, its 
geometry enhances its importance. Separating the kx~ 
and fcz-dependence arising from the terms Ji = 2Ji — J[, 
J", J2, J2, and J2 is in principle a complex task, and tun- 
ing of these parameters would permit a very accurate fit 
to the measured dispersion. For the purposes of the cur- 
rent, unrefined fit we make use of the remarkable observa- 
tion that the magnon energy in direction [1,0,-2] (line AC 
in Fig.|21for TlCuCls; a similar observation was made for 
KCUCI3) is small and disperses only weakly and mono- 
tonically. In a tight-binding model this may be achieved 
only if the [k^^kz) dispersion is largely determined by 
the term Jl^. As in Refs. andU^ we find that a good 
fit can be achieved to all of the dispersion features using 
only the parameters Ji, J2 and J3, with contributions 
from the other parameters being negligible. From this 



procedure it is not possible to distinguish between the Ji 
and J[ contributions to Ji , but the result that this is neg- 
ative is robust. This implies that the Ji exchange path 
is closer to the FM coupling regime achieved near Cu- 
Cl-Cu bond angles of 90° than is the J[ path [Fig. 1(a)]. 
Although it is likely that in this case both values are 
individually rather small, the bond-operator framework 
provides a simple origin for the near-cancelation of two 
AF interactions to give a small net coupling despite the 
short bond lengths. 

In such a scheme, the exceptionally large required 
to fit the dispersion of TlCuCls (Table I) has no imme- 
diate justification from qualitative arguments based on 
superexchange paths. However, it is also by no means 
implausible given the tilted dimer orientation within the 
unit cell (Fig. 2) , which results in all of the bonds on this 
path being rather straight. However, because this cou- 
pling provides the only cos 2kx modulation, and the only 
cos kz modulation, its value is determined very accurately 
by its contribution to the dips away from the dispersion 
maxima towards the zone boundaries in both k^- and 
/c^-directions, primarily in competition with J3. For the 
strongly dimerized, 3d systems represented by both sets 
of parameters in Table I, these effects may only be due 
to competition between interactions of different periodic- 
ities, and there is no possibility of such dispersion forms 
being provided purely by many-body effects of the type 
expected in weakly coupled spin chains»^ Such competi- 
tion is ensured for the assumed dimer coupling geometry 
by the presence of the pathways within J3. 

As with Ji, the zero-field data alone does not per- 
mit the possibility of fitting both contributions J3 and 
Jg to the effective coupling J3. However, because the 
minimum of the magnon bands is located at (0, 0, 2ti) in 
reciprocal space, one may expect (on the assumption that 
all couplings are AF) that J[ > 2Ji, whence the nega- 
tive signs for the effective parameters Ji in Table I, and 
that J3 > J3. One attempt has been mad©^ to extract 
these parameters for TlCuCls using a cluster expansion 
theory)^ which returns the intersite couplings (valid for a 
formulation different from that of Sec. II) Ji = 0.34meV, 
J[ = 1.70meV, J3 = 0.91meV, and = -0.57meV. 
However, for the present purposes, namely a focus on 
the ordered phases of the system under applied field and 
pressure, we choose to avoid the large uncertainties in- 
herent in these estimates and proceed simply with the 
effective quantities Ji and J3, determined at zero field, 
in place of iJi and J3 in Eq. ||SJ|, with Tii and set 
to zero. This approximation is in fact insufficient for 
a full treatment of the ordered phase, where the coeffi- 
cients of the terms TLtt have the same sign, and should 
no longer be necessary when reliable high-field measure- 
ments become available. The resulting small deviations 
from quantitative accuracy which are found in reproduc- 
ing the high-field magnetization (Sec. HIE) imply that J[ 
and Jg are, respectively, significantly larger than Ji and 
Jg. We have thus extracted only four exchange constants 
(Table I) as the minimal model required to describe 
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TlCuCls and KCuCls, and have shown in Fig.Othat this 
set provides a remarkably good account of the magnon 
dispersion relations in both systems. In the following 
sections we therefore employ this model with the param- 
eters, ( J, Ji, J2 , J3), and within the Holstein-Primakoff 
approximation (Sec. IID), to study field- and pressure- 
induced magnetic ordering in TlCuCla and KCuCla. 

We comment here that the XCuCla structure is not 
magnetically frustrated: when following the dominant 
interactions no triangles of spins are found. The mod- 
ulation of the dispersion relations arises rather from the 
simple fact that components of different periodicities may 
contribute as a consequence of the complex unit cell and 
coupling geometry. The lack of frustration is a very 
important feature for both the appearance of staggered 
magnetic order in an applied field and for the success of 
the bond-operator description. 

With regard to a comparison of the parameters in Ta- 
ble I with the results of other approaches, these are in- 
deed rather similar to the values given in Refs. 4 [differ- 
ences are due to the fact that the experimental fitting for- 
mula is linear, whereas Eq. H16() has a square-root form] 
and The bond-operator formulation has the advan- 
tage of providing a direct justification for AF interdimer 
bonds J2 and J3. The opposite sign of the term found 
by the alternative technique is rather harder to justify on 
the grounds of magnitude and orbital overlap geometry. 
The differences between the results of Table I and those 
shown in Table I of Ref. [s^ are due solely to the con- 
ventions chosen for the bond-operator analysis. That 
the magnon excitations are so well described by such a 
simple model indicates a posteriori that the qualitative 
assumptions made in Sec. IIA concerning the relevant ex- 
change paths are appropriate. The increase in interdimer 
couphng strengths observed from KCuCla to TlCuCls, as 
the anion size increases, may suggest that the role of the 
anion in determining the superexchange parameters may 
depend more on its charge distribution (active participa- 
tion in superexchange) than merely on its spacing effect. 

In summary, TlCuCla and KCuCla are well described 
by simple models of dimerized Cu spins coupled by rel- 
atively weak superexchange interactions in all three spa- 
tial dimensions. This model may be used not only for 
magnon dispersion relations measured by INS, but also 
for the intensities of the elementary excitations, which are 
dominated by the dimer form factor. With these zero- 
field results we turn now to a treatment of the excitations 
in magnetic fields H > He sufficient to close the spin gap 
in TICUCI3. In Sec. IV we will return to these consid- 
erations to justify our treatment of the superexchange 
integrals under applied hydrostatic pressure. 



III. FIELD-INDUCED ORDER 

In the previous section we have established a mini- 
mal model within the bond-operator framework which 
describes the ground state and magnetic excitations of 
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TlCuCls and KCuCls at zero magnetic field. In this sec- 
tion we will apply the same formalism to consider the 
evolution of the system under an applied field, with par- 
ticular focus on the excitations of TlCuCla in the exper- 
imentally accessible intermediate-field regime H > He, 
where the lowest-lying magnon mode is gapless. 

The bond-operator formulation for ground and excited 
states in magnetically ordered phases was first explored 
for systems where the order arises as a result of increasing 
the superexchange interactions between low-dimensional 
subsystems with a spin gap.-— In this case, all three 
magnon modes in the disordered phase remain degener- 
ate as the coupling rises and the spin gap is closed, until 
at the quantum critical point (QCP) the three magnons 
become massless at the band minimum (magnetic Bragg 
point). On further increase of the coupling into the or- 
dered phase, two of the magnons remain massless, corre- 
sponding to spin waves, while the third acquires a mass 
which increases on moving away from the QCP. This last 
excitation has the interpretation of a longitudinal mode, 
and corresponds to fluctuations in the amplitude of the 
order parameter, which are soft close to the QCP. This 
situation will be discussed further in Sec. IV in the con- 
text of the pressure-induced QPT in TlCuCla. 

The situation in a magnetic field is rather different, 
and has been considered using bond operators for the 
isolated spin ladder;^ for a pair of coupled planeS)^ and 
for TlCuCla in Ref. ^Ed. A simple understanding may be 
obtained by considering a single, isolated dimer with AF 
coupling J: at zero field the ground state is a singlet, and 
the three-fold degenerate triplet excitations at energy 
A = J separate in an applied field due to the Zeeman in- 
teraction, as shown schematically in Fig.QJa). At a criti- 
cal field, given hy he = A, the energy of the lowest triplet 
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is reduced to zero, and the crossing of levels changes 
the ground state from singlet to triplet. Above He, the 
ground state is of pure triplet nature and the lowest-lying 
excitation is a singlet. For a set of dimers with mutual 
AF interactions [Fig. 2fb)], the zero-field triplet excita- 
tion modes are dispersive in reciprocal space, and the 
critical field He is defined by the value required to close 
the zero-field gap, at which point a single mode becomes 
massless. 

The ground state of the ordered phase has the two 
essential properties that is it possesses a field-induced 
staggered moment and has a massless excitation corre- 
sponding to the transverse fluctuations of this moment. 
The appropriate description of this phase is a linear com- 
bination with the singlet of the lowest- and highest- 
lying triplet states,^" which corresponds to the choice of 
a mean field for the staggered moment, whence the trans- 
formation in Holstein-Primakoff formulation returns the 
excitations of the ordered phase. We note here that the 
conventional bond-operator treatment (Sec. IIC) may 
lose consistency in the description of these excitations be- 
cause higher-order interactions between triplets are con- 
tained in an uncontrolled manner. The coefficients of the 
ground-state singlet-triplet admixture change continu- 
ously with the applied field until the upper critical field 
Hs, where all spins are aligned. Beyond Hg, the ground 
state has pure triplet character (Fig. ^ and three ex- 
cited modes of the strong-field regime are all massive, as 
expected for a field-aligned antiferromagnet. 



A. Disordered phase 

Henceforth we adopt the notation Ji for Ji, J2 for J2 , 
and J3 for J3. As in Sec. IIC, the dispersion relations for 
the triplet magnons of the disordered phase are given by 

Eka = VJ^ + 2Jek " ah, (18) 

Efc ~ — i Ji COS fcj, — i J2 COS (2fc2, -I- fc^) 

+J3 COS {kx + ^kz) COS ^ky, (19) 

and are shown in Fig. |31for the critical field he, although 
a qualitatively identical picture is obtained for all fields 
h < he- We emphasize that the lowest mode retains its 
quadratic shape at the band minimum (c/. Sec. IV) even 
at he, while beyond this, as we will demonstrate below, 
an additional linear component arises. 

The nature of the three magnon modes in the disor- 
dered phase may be understood as follows. The wave 
function of the lowest triplet mode, E^._^, can be approx- 
imated by a linear combination of singlets and triplets on 
each dimer bond, 

^^ ~ ^"(1 Ti) - I Ti)) - «e'('=-"'-^fc+*) I TT), (20) 

where u is a coefficient of order unity and v is small and 
real. The expectation values of the spin operator com- 
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FIG. 5: Magnon dispersion for TICuCls aX H ^ (= 5.6T), 
using g — 2.16 for the ^-factor of TlCuCla. 

poncnts at a given dimer i in the state (|20|l are 
{Si,x)i = -{Sr,x)i ^ ^uvcos{k ■ ri-Ej^j^t), 

{Sl,y)i = -{Sr,y)t ^UV Sm{k ■ Ef^^t) , (21) 

{Sl,z)i = {Sr,z)i ^ \v^, 

from which it is evident that the mode possesses a very 
small, uniform, static magnetic moment parallel to the 
field, and thus gains Zeeman energy. Perpendicular to 
the field there is also a finite magnetic moment, which 
fluctuates in both space and time with characteristic 
wavevector k and energy i?fc_|_, and is staggered (spins 
I and r oppositely aligned) because of the AF intradimer 
coupling J. Similar considerations for the mode E^._ 
reveal that the uniform magnetic moment is directed an- 
tiparallel to the field, leading to a higher Zeeman energy. 
Finally, the E^.^ mode has no magnetization perpendicu- 
lar to the field, although it does possess a moment paral- 
lel to the field with modulation wavevector fc, and so its 
energy does not change with field. For H < He, the en- 
ergies of these modes are positive, meaning that the spin 
structure associated with the magnons is very weak and 
fluctuating. We stress that the sum of the magnetiza- 
tion components remains zero in the disordered phases, 
which can be considered as a consequence of quantum 
fluctuations which act to restore the 0(3) symmetry of 
the field-free system. 

We note again that an increase in the applied field leads 
to a shift in mode energies without changing the shape of 
their zero-field dispersion, until the lowest mode {a = -\-) 
becomes soft at Q = (0, 0, 2tt). This determines both the 
critical field, 

he - gimHe = v/J2 - J( Ji + J2 + 2J3), (22) 

and the spin structure of the lowest-lying mode, given 
by fe = Q, which corresponds to the static spin con- 
figuration in the field-induced ordered phase. One may 
conclude that a TlCuCla system of classical spins would 
favor AF order to minimize simultaneously all intersite 
interaction terms. By contrast, for the quantum system 
the strong AF intradimer interaction, J, favors dimer 
singlet formation, and the AF spin structure cannot be 
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the ground state for H < He- However, the incipient, 
long-ranged AF spin order is reflected in the structure 
of the lowest-lying magnon excitation, which is minimal 
aX k — Q. As the magnetic field is increased, the lowest 
mode is lowered progressively until i/c, where dimer sin- 
glet formation is no longer sufficient to cause an ordered 
phase, the symmetry of the system is lowered to 0(2), 
and the onset of real uniform and staggered magnetiza- 
tion components is observed. For the magnon excitations 
in the ordered regime, the instability is signaled by the 
lowest mode becoming negative within the formulation of 
this subsection. We turn now to a consistent description 
of all of these phenomena in the ordered phase {H > He). 



B. Ordered phase {He < H < Hs) 

Magnon excitations 

In the intermediate-field regime, the ground state of 
each dimer can be considered as a partially polarized FM 
configuration. It is important to note that there is no ex- 
plicit AF component in the dimer ground state, and that 
the ordering emerges only from closing of the gap to the 
lowest magnon branch aX k — Q. In the bond-operator 
formulation, this ordered ground state is represented^2iS 
in the basis of Eq. (0J by finite expectation values s , t+ , 
and i_ of these singlet and triplet operators. The com- 
ponent of the highest-lying triplet mode in the ground- 
state condensate may appear counterintuitive, and was 
neglected in a number of approximate treatments i^SiiS 
However, the presence of terms of the form SjSi 

makes clear 



Holstein-Primakoff representation to 



in the transformed Hamiltonian of Eq. Q 
that a finite component of this state is required in the 
consistent condensate. 

Following the treatment of the preceding sections, 
the linear combination of operators appropriate for the 
ground state mode is taken as the condensate, and the or- 
thogonal mode containing a singlet component describes 
its gapless phase fluctuations. For a full description of 
the ordered regime we perform the transformation^ 



= tif)^ (23) 

h- = fU- - gU+, 

in which the fc-independent coefficients u, v, /, and g 
arise from successive unitary transformations. The con- 
ditions u'^ + v'^ = 1 and P+g^ — 1 dictate that u — cos 9, 
V = sin0, / — cos(f), and g = sin0, with 9 and (p the two 
independent variables required to specify the condensate. 

As above, the operator is treated as a uniformly 
condensed mean-field parameter, 

a, a. (24) 

Because the two unitary transformations conserve parti- 
cle number, the global constraint may be reduced in the 



\ 



- E 



(25) 



We emphasize again the participation of the highest 
triplet mode (t_) in the condensate, as a consequence 
of the triplet pair-creation processes in Hst © ■ The lin- 
ear combination of singlet and triplets in the condensate 
a yields a staggered magnetization perpendicular to the 
field, whose order is given by the wave vector Q. 

Using the transformed operators H23|) . and with re- 
placement of the condensate parameter using the global 
constraint of Eq. 125|) , the Hamiltonian © is reexpressed 
up quadratic order in the b operators in the form 



H^O{b^) + 0{b^) + 0{b^), 



(26) 



in which the classical part at zeroth order is given by 

O(feO) = N{{e+f + e^g^y (27) 
-eQ[-u\,\l + 2fg) + lvXf-gY]}, 

where 



eg 



J — ah, (a : 



±) 



(28) 
(29) 



The terms at linear order may be written as 



e-5 

.2/ i-2 



0{b') = w«(46Q+ + H.c.){e+/2 

+e^Q[{u'^v'){l + 2fg)~v'if' 
+uv{albQ_ +II.C.) (30) 
X [(-£+ + e_)/.g + eQ(«2 + 2fgv^){f - g% 



and the requirement that this sum be zero fixes the pa- 
rameter choice {9,(f)), which also minimizes 0(&°). The 
equations determining the two unitary transformations 
are then given by 

e+f + e-g' = ~eQ[iu'-v')il + 2fg)-v'if~gYl 
(-e+ + e_)/5 = -eQiu^ + 2fgv^)if~g^). (31) 

The field dependence of the coefficients u^, w^, and g^ 
is shown in Fig.^lfor the parameters of both KCuCls and 
TlCuCls. Below the critical field {H < He), the ground 
state is a pure singlet condensate specified by u = 1, 
V ~ 0, f — 1, and 5 = 0, i.e. by the angles {9, </>) = (0,0). 
Precisely at the critical field, where u = 1 and v = 
{9 — 0), g takes on a finite value while / < 1. In the 
intermediate regime, 9 then changes continuously from 
to 7t/2. Except for the parameters of TlCuCla and for 
fields close to He, (f> remains small, reflecting the fact 
that mixing of the upper and lower triplet modes is not 
generally strong. Finally, above the saturation field {H > 
Hg), the ground state becomes a pure condensate of the 
lowest triplet, which is described byu = 0,w = l,/ = l. 
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FIG. 6: Magnetic-field dependence of tlie coefficients , , 
f, and g\ for (a) KCuCla and (b) TlCuCla. 



and g — 0, or by {0,(p) = (7r/2,0). The critical field is 
determined exactly as 



he = gtiBHc = ^.P- J{Ji + J2 + 2J3), 



(32) 



which coincides with the soft -mode condition (|22|l in the 
low-field regime. For TlCuCls (KCuCla), the parameters 
of Table I give a critical field = 5.6T (19.6T), con- 
sistent with the measured valuejSi The saturation field is 
determined by the condition 



gfiBHs = J + J1 + J2 + 2 J3 



(33) 



required to overcome all of the AF bonds to a single spin. 

The spin excitations in the ordered phase are described 
by the quadratic term in Eq. (|26l) . 



oib^) = no + n±, 



where 



k 



(34) 



(35) 



and 



k 

+ (efc±^fc+&fc-+H.c.) 

+ i(A^_6fc_5_fc_+H.c.) 
+ (A^^6fc+6_fc- +H.C.) 



(36) 



For continuity of presentation, the expressions for e^j and 
Afcj are given in the Appendix [Eqs. (|A.ip and l|A.3p ]. We 
note that on setting (0, (f) = (0, 0) the Hamiltonian [|23 
and (|34|) ] reduces to that of the disordered phase given 
in Eq. (fT5|l with 's — 1 and = — f J. 

The 0{b^) terms in Eq. H34|) are diagonalized by two 
separate Bogoliubov transformations which yield the en- 
ergies of the collective modes of the condensate. Because 
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FIG. 7: Dispersion of the three magnon modes in TlCuCls: 
(a) at H = 14T > (cf. Ref.*!); (b) at the saturation field 



the triplet b^a does not mix with the other two modes, 
Tio is diagonalized easily to give 



-E'fco 



'fco 



A2 



(37) 



which describes the middle magnon branch in the ordered 
phase. TL± is diagonalized by the transformation 

«fc = + n^i:-bk_ + i;^Vfe-fc+ + (38) 

where m = 1, 2 denotes the two modes corresponding to 
each value of fc, whence 



mt n 



k "fc 



(39) 



and are then the eigenvalue and eigenvector for 
the mode of branch to, and are related by the equation 




f^k+\ 



VkJ 



(40) 

This equation has four eigenvalues, of which two are posi- 
tive, and we define the smaller of these as E'^, the disper- 
sion relation of the lowest magnon mode, while the larger 
is defined as E^^ for the highest mode. By the nature 
of the problem, the other two eigenvalues are —E^ and 
—EJ^. In the disordered phase, efc-t = = A^^ = Afc_, 
and the equation H40|l reduces to a 2 x 2 matrix form. 
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FIG. 8: Field dependence of the energy gaps for the three 
magnon modes in TlCuCls: (a) points are data from the INS 
experiment of Ref. ;iS, solid lines are the theoretical fit using 
the zero-field parameters of Table I; (b) field-dependence at 
high fields, using g = 2.16 as the g-factor of TlCuClai^ 



In the ordered phase above He, the mode Ef._^. is al- 
ways massless at the band minimum, Q, and develops a 
linear dependence on k (in all three directions in recip- 
rocal space) in the vicinity of this point [Fig. IJfa)]. This 
branch is the Goldstone mode, which arises from the fact 
that rotation of the staggered moment about the axis of 
broken symmetry (the field axis) does not change the en- 
ergy of the system. It is the staggered magnetic order in- 
duced perpendicular to the applied magnetic field which 
breaks rotational symmetry around the field axis. Math- 
ematically, rotations of the induced moment are realized 



by changing the phase of / and g in H23|l (/ 



'1-X 



/, 



g — !• e^^g, where x is the rotation angle), and the invari- 
ance of the expressions under this transformation deter- 
mines the gapless nature of the Goldstone mode. The 
remaining modes in Fig. Cfa) show that for fields close 
to He the overall shape of the dispersion relations is not 
strongly affected by the presence of the field while the 
field remains small on the scale of the mode energies Ef,. 

The field-dependence of magnon dispersion relations 
may be characterized up to high fields using the expres- 
sions for obtained for E/^q, and with the zero- 
field parameters of Table I. Only in the disordered phase 
{H < He) does the magnon dispersion not change its 
shape, with a simple linear splitting of the three branches 
developing due the Zeeman interaction until the energy 
gap of the lowest mode is reduced to zero at H = He- 
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FIG. 9: (a) Normalized magnetization curves for KCuCls. 
Open circles are experimental data for Mz measured at T = 
1.3K»itt The theoretical fits are made using a = 2.29 as the 
g-factor for KCuCls.^" (b) Normalized magnetization curves 
for TlCuCla . Open circles are experimental data for Mz mea- 
sured at T = 1.3K(- and closed circles are data for M^y mea- 
sured at T = 0.2Ki— The theoretical fits are made using the 
same g-factor for TlCuCla as in Fig. |5] 



Above the critical field, the gap for E^.^, the Goldstone 
mode of the ordered phase, remains zero, as a conse- 
quence of which the gaps of the higher modes, EgQ and 
Eg_, show an abrupt increase in slope at He [Fig.|HIa), 
cf. Fig. 01]. The results of the bond-operator theory 
(solid lines) show quantitative agreement with the ex- 
perimental measurements obtained by INS.^ 

Fig. IHIb) shows the predicted field dependence of the 
gaps for high fields: while Eg^ has monotonically increas- 
ing behavior all the way to and beyond Hs, the gap EgQ 
of the middle branch first increases with field, but then 
decreases until Hg. This behavior is related to the posi- 
tion in the Brillouin zone of dispersion minimum, which 
is determined largely by the e^Q term in the Appendix. 
As the magnetic field is raised, the sign of the dispersive 
component changes on passing through the field at 
which = v'^, which for the parameters of TlCuCls cor- 
responds to 55T [Fig. Elb)]. Below this value, where 
> w^, the minimum is at Q, while for fields close to 
55T {u^ ~ w^) the dispersion becomes flat. Above this 
field (it^ < u^) the minimum moves from Q to a point in 
the middle of the zone, as shown for the branch E^q in 
Fig. Hb) , and Ego then decreases with field until satura- 
tion. Finally, at H = Hg, the k dependence of the lowest 
mode, becomes quadratic again in the vicinity of 

Q [Fig. Ctb)] as a gap opens to the lowest mode of the 
fully saturated ferromagnet [Fig.|SIb)]. For H > Hg (see 
below), the magnon modes decouple, the condensate has 
pure triplet character and the lowest excitation mode is 
a pure singlet whose gap grows linearly with the field, 
as do EgQ and (with a slope greater by a factor of two) 
Eg^. The evolution with field of the shape of the dis- 
persion relations through the ordered phase is illustrated 
by Fig. Ul where it is clear that the lower branch changes 
rather little, the upper branch becomes progressively nar- 
rower until it is completely non-dispersive at Hg , and the 
middle branch is inverted as discussed above. 



12 



Magnetization 

The theoretical framework also gives direct access to 
the average magnetization of the dimers in the ordered- 
phase condensate. We denote by the magnetization 
paraUel to the field, whose normalized value is 

^'i^^-^Y.^A+U+-t\-U-). (41) 

i 

and by Mxy the perpendicular (staggered) magnetization 

i 
i 

where i labels the dimers and (...) signifies an expecta- 
tion value calculated using the operators in Eq. (|^ . 
The mean-field magnetizations, normalized to their val- 
ues at saturation are given by 

= v^if-g"). (43) 
M^y = V2uv{f + g). (44) 

These expressions correspond to contributions, and 
are supplemented by terms of the form a/^a'^, which 
correspond to quantum corrections due to magnon ex- 
citations. Both contributions are shown in Figs. Efa) 
andinib), but the quantum corrections are found to be 
small, presumably as a result of the three-dimensionality 
of the KCuCla and TlCuCla systems. In both cases the 
square of the staggered moment (M^y) shows a linear 
field dependence close to He and to Hs, indicating that 
M^y \/ H — He as expected from a mean-field descrip- 
tion. 

For KCuCls, the magnetization parallel to the field 
(Mz) is almost perfectly linear in H [Fig.|^a)], in good 
agreement with experiment.^'' If one neglects the upper 
triplet mode (t-) in the condensate specified by Eq. (|23|l . 

would be completely linear in H, and because all in- 
terdimer interactions are small (Table I) , KCuCla is close 
to this limit. By contrast, the interladder interaction J2 
is considerably larger for TlCuCls, where mixing into the 
condensate of the i_ mode, mediated primarily by inter- 
actions of the type t^iLss, is significant. This point can 
be seen clearly by comparing Figs.E^a) and|SJb), where 
the value of is considerably larger for TlCuCls than for 
KCuCla in the vicinity of He- Physically, the presence of 
the highest triplet in the condensate costs Zeeman energy, 
which leads to a suppression of v, and consequently 
is reduced near He for TlCuCls. This causes a significant 
deviation (due in fact to higher powers in H — He) from 
purely linear field dependence, as illustrated in Fig.|^b), 
which is again in very good agreement with the observed 
form.^ This result emphasizes the importance of includ- 
ing the t- triplet in a consistent description of the ground 
state, which is required to account for the difference be- 
tween the magnetization curves of TlCuCls and KCuCla. 



C. Saturated phase {H > Hs) 

In the high-field regime the spins are fully polar- 
ized. In the notation of Sec. IIIB, the saturated phase 
is described by {0,4)) = (7r/2,0) in the Hamiltonian of 
Eq. which specifies that the ground state becomes 

a pure condensate of the lowest triplet (t+ in the notation 
of Sec. IIIA). This phase is characterized by a complete 
decoupling of the triplet modes, by a finite energy gap 
to the lowest mode, which corresponds to an excitation 
from the triplet ground state to a singlet, and by the 
fact that the highest triplet mode becomes dispersionless 
[Fig. [7{b)]. In this regime the interaction term Htt be- 
tween the lowest and highest triplet branches in Eq. Hll() 
becomes an on-site potential term for i_ magnons, and 
does not assist their propagation, so there is no first- 
order hopping process for this branch. A possible second- 
order process exists for propagation of the upper triplet 
by creation of two singlets in a virtual process due to the 
term sjs]ij+ti_ in Hst- The quantitative effect of this 
term on the dispersion is expected to be very small, 
but the process may also have a qualitative role in me- 
diating the decay of t_ magnons into two singlet modes 
(bearing in mind the triplet nature of the condensate). 
With energy and momentum conservation specified ac- 
cording to 

Ek- = J2^Ek+g,+ + Ek-q.,+), (45) 

this may lead to a finite lifetime for the highest triplet 
excitations. We note also that this process is possible 
for fields He < H < Hs, and provides one reason for 
which the upper magnon mode may become broad at 
high magnetic fields. 



IV. PRESSURE-INDUCED ORDER 

In Sec. Ill we have discussed field-induced magnetic 
order in a quantum spin system, which can be under- 
stood qualitatively as the consequence of the lowering of 
one of the triplet modes by the applied magnetic field to 
the point where it becomes soft, and the ordered phases 
is favored. Another way in which the magnon modes may 
soften, causing the onset of magnetic order, is by increas- 
ing the interdimer interactions, broadening the magnon 
band. The application of pressure is expected to have just 
this effect, and for the small-gap system TlCuCls this 
type of pressure-induced AF order was found recently by 
Oosawa et al}"^ in neutron diffraction experiments con- 
ducted under a fixed hydrostatic pressure P = 1.48GPa. 

As noted in Sec. Ill, this type of "pressure-induced" 
magnetic QPT provided the first application of the bond- 
operator formulation to ordered phases »2i In the disor- 
dered phase the three magnon modes are degenerate at 
H — for all couplings, including at the QCP where 
they are gapless at the wave vector Q characterizing the 
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incipient magnetic order. The excitations of the ordered 
AF consist of two spin waves, which are the conventional 
phase modes (Goldstone modes) of the ordered moment, 
and a third which corresponds to its amphtude fluctua- 
tions and is generaUy not relevant for the low-energy dy- 
namics, but may become low-lying close to the QCPiSl 
The ordered phase is described in the bond-operator ap- 
proach by a condensation of one of the triplet operators, 
conventionally taken as t| in Eq. iQJ , and in Ref . this 
was treated as a separate condensate in order to obtain 
an effective description of the ordered phase close to the 
QCP. 

A complete account of the changes to the superex- 
change interactions caused by the application of a hy- 
drostatic pressure would require a detailed structural in- 
vestigation of the alterations to bond lengths and angles. 
Such a study lies beyond the scope of the current analysis, 
the aim of which is to characterize the changes in phys- 
ical properties through the QPT. We thus employ qual- 
itative arguments based on the considerations of Sec. II 
for the probable evolution of the coupling constants in 
the XCuCla structure, and use this as the foundation for 
a plausible model. 

The structural unit composed of CuClg octahedra is 
not expected to undergo significant changes except at 
rather high pressures, and so we assume here that the 
dimer coupling J is unaltered. The leading effects of 
hydrostatic pressure are expected to be a small longitu- 
dinal compression of the ladder units, and, due to the 
absence of strong covalent bonds in the transverse direc- 
tions, a somewhat larger reduction in their separation. 
These changes generally involve small reductions in bond 
lengths and more significant distortions of the bond an- 
gles. Considering next the effective interdimer coupling 
Ji = 2Ji — J[ (in the notation of Sec. II), any small in- 
creases in ladder leg coupling Ji are likely to be offset to 
a considerable extent by similar changes in the diagonal 
coupling J[ , so we will also neglect the pressure-induced 
alteration of this parameter in the following simplified 
treatment. The leading pressure effects are thus expected 
to involve the alteration of the parameters J2 and J3, 
both of which are determined by Cu-Cl-Cl-Cu path- 
ways, although the latter may also have contributions 
from the X+ ion. We assume (based on the experimental 
observation of the QPT) that J2 and J3 increase with the 
pressure, P. For the functional form of this increase, the 
exact dependence of superexchange interactions on bond 
lengths and angles remains a topic active research, and 
in the absence of structural data under pressure we will 
attempt to fit this form with minimal assumptions. We 
take the interladder interactions to change with applied 
pressure according to 



J2(P) = p(P)J2, 
UP) = p{P).h, 



(46) 



where p{P) is a dimensionless function of the real pres- 
sure P and the values of J2 {J 2) and J3 ( J3) at p = 1 
are listed in Table I. We begin by analyzing the situation 
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FIG. 10: Threefold degenerate triplet magnon dispersion at 
zero field. Dashed and solid lines correspond respectively to 
atmospheric {p = 1) and critical {p = Pc — 1.018) pressures. 



where the interladder interactions are simply increased 
by a factor p, and then discuss the experimental deter- 
mination of the hmction p{P). 



A. No magnetic field 

With a view to the simultaneous application of field 
and pressure (below), we discuss the pressure-induced 
QPT using the triplet operators t+, to, and t- which 
were introduced for the field-induced QPT, rather than 
in the "natural" basis of operators t^, ty, and In 
the absence of a magnetic field, the magnon dispersion 
in the disordered phase is given by Eq. H16|) with h — 0, 
and the single branch is triply degenerate. The disper- 
sion, shown at atmospheric pressure by the dashed line 
in Fig.^(c/. Fig. El, is given by 



(47) 



where e^, is defined in (|19|l but with J2 and J3 now 
given by Eq. 146|) . As pressure is increased, the maxi- 
mal positive and negative values of are enhanced, the 
bandwidth increases and the dispersion minimum (spin 
gap) decreases. It is clear from Fig. ^3 that the low- 
energy part of the magnon spectrum is particularly sen- 
sitive to the pressure, whereas the higher-energy parts 
are scarcely affected. This result is a simple consequence 
of the square-root form of Eq. (|47|l , which amplifies the 
changes in E/^ where the two terms are close to mutual 
cancellation, i.e. in the vicinity of the band minimum Q. 

The magnon modes become soft, Eq = 0, at a critical 
pressure, Pc, where the enhanced interladder interactions 
favor a long-ranged AF spin order over a gapped dimer 
phase. The dimensionless pressure p{P) at the critical 
pressure for the parameters of TlCuCls is 



Pc 



J - Ji 
J2 + 2 J3 



1.018, 



(48) 



from which one notes immediately that the critical pres- 
sure represents only a very small fractional increase in 
P- 
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FIG. 11: Magnon dispersion relations in pressure-induced 
AF phase at zero field, showing two degenerate phase modes 
(solid line) and one amplitude mode (dashed line). 
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FIG. 12: Pressure dependence of spin gaps for phase and 
amplitude modes. 



Expansion of in Eq. (|47|l around Q, CQ+^fc 
(Jefc, yields a dispersion 



Q+Sk 



2J5ek. 



(49) 



Because contains cosine terms in k which are min- 
imal at Q, 6ei^ is quadratic in Sk, and at the criti- 
cal pressure Ef^ has a linear fc~dependence around Q 
[Fig. ^]. This stands in explicit contrast to the case of 
field-induced order, where the single low-lying magnon 
retains a quadratic dependence on k at the critical field 
where it softens. 

Increasing the pressure beyond pc induces the AF or- 
dered phase, which may be described by the same oper- 
ator transformation as in Sec. III. In the absence of an 
applied magnetic field, the coupled triplets t+ and i_ are 
equivalent and / = g = 1/V2 7r/4) in Eq. (j2SJ- 

We note here that this combination in fact restores the 
eigenvectors and ty in Eq. (QJ , referred to above as the 
"natural basis" for the field-free problem. Using the b op- 
erators defined in Eq. H23|) , the Hamiltonian at quadratic 
order takes the symmetrical form 



Ho + 'H± — TLa 



(50) 



with 



= 5I[efca&L^fc" + U^icPko.b-ko. + H.c.)]. (51) 
k 

The matrix elements are given by 

efco = efc- = J"^"^ - 4eQU^w^ + [u^ - v'^)ek, 
Afc+ = [u'-v'fek, (52) 



where 



2 1 J 
u = - ' 



1 J 



2 Meg)' 



2 (-4eQ) 



(53) 



and 

eQ = -i[Ji+p(^2 + 2J3)]. (54) 
The magnon mode dispersion relations are then given by 



ka 



ka ka ' 



(55) 



and are illustrated in Fig. ^2 With this choice of coordi- 
nates the staggered magnetic moment is developed along 
the X-axis in spin space, as the transformed magnon op- 
erator bff.^ is precisely tx of Eq. . It is clear that the 
other two modes remain degenerate, and that E/^q and 
are gapless at Q with linear dispersion around this 
point. These are the spin waves of the ordered magnet, 
or Goldstone modes; their magnetic moment fluctuates in 
the plane perpendicular to the induced moment, demon- 
strating their role as the transverse phase modes. The re- 
maining mode, has a finite gap and corresponds to 
fluctuations in the direction of the induced moment, con- 
firming that this is the longitudinal, or amplitude, mode 
of the pressure-induced AF phase. This is the splitting 
of the threefold degenerate triplet modes expected from 
previous studies of coupling-induced AF order in gapped 



spin systems 
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The excitation gap of the amplitude 



mode may be written as 



E, 



2uvy/2 



(56) 
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and its evolution with applied pressure is shown in 
Fig.im 

We turn next to a discussion of the pressure de- 
pendence of the interladder superexchange parameters, 
i.e. to a determination oi p{P). The staggered moment 
grows with increasing pressure in the ordered phase as 
shown in Fig. 1131 and close to the critical pressure it has 
the mean-field exponent M,j.y oc {p — PcY^^ . Unlike the 
field-induced order, there is no uniform magnetization 
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P 

FIG. 13: Pressure dependence of induced staggered moment, 
Mxy, in TlCuCla at zero field. The normalization is to a 
saturated moment of 1. 



thus be performed at lower pressures and with larger 
samples. TlCuCla therefore offers an ideal system for 
observation of the amplitude mode, and for the measure- 
ment of its dispersion over a significant portion of the 
Brillouin zone. These measurements could be used to 
test the hypothesis that the amplitude mode is broad, 
and potentially ill-defined, due to the possibility of de- 
cay processes into pairs of spin waveSfS4 The gap of the 
amplitude mode may also be measured from the loca- 
tion of additional contributions to susceptibility, nuclear 
magnetic resonance (NMR), and electron spin resonance 
(ESR) signals. 



at zero field, which is reflected in the theoretical descrip- 
tion by the relation f — g. The critical pressure has been 
found to be approximately Pc — 2kbarii& By expanding 
the function at the critical pressure, a crude determina- 
tion of p{P) may be made by assuming the form 

p{P)-p, = P^{P~P,)+P2{P~Pcf. (57) 
where /3i and [32 are constants to be determined from 

Vo-Pc = MPo ~ Pc) + HPo - Pc)^ (58) 
Pi-Pc = Pi{Pi-Pc)+l32{Pi-Pcf- (59) 

The three known data points are i) atmospheric pres- 
sure P ~ Pq, where pq — I, ii) the critical pressure 
P = Pc = 2kbar, where p = pc = 1.018, and Hi) 
P = Pi = 1.48GPa (14.8kbar), where the staggered 
magnetization M^y has attained 60% of its saturation 
value, which corresponds to pi = 1.4 (Fig. I13|l . With 
these values the prefactors of the lowest orders in the ex- 
pansion are f3i = 0.012kbar"^ and P2 = 0.0014kbar"^. 
The relative decrease in magnitude of these coefficients 
suggests that higher orders may not be significant, al- 
though a detailed experimental determination is awaited. 
In the vicinity of the critical pressure, the finite linear 
term in the pressure dependence of p{P) leads to an ex- 
ponent M^y (X (P — Pc)^^^ . The most important qualita- 
tive observation arising from these considerations is that 
obtaining a critical pressure of Pc = 2kbar is straight- 
forward in the laboratory, which in practice implies that 
experiments probing the properties of the TlCuCla sys- 
tem on both sides of the pressure-induced QPT could be 
performed on samples of considerable volume. 

The pressure dependence of the spin gap of the am- 
plitude mode (|56|l . shown in Fig.^lis similar to that of 
the induced staggered moment fFig. ll3() . as both depend 
in a similar manner on the singlet-triplet mixing coeffi- 
cient 6. Under a pressure of 1.48GPa, where the neutron 
diffraction experiment revealing magnetic order was per- 
formed, we estimate the gap of the amplitude mode to 
be 5meV. This is well above the value required to dis- 
tinguish the amplitude mode from the spin waves and 
magnetic Bragg peaks in INS experiments, which may 



We conclude this section by noting the two qualita- 
tive differences between the pressure-induced QPT and 
the field-induced situation discussed in Sec. III. The first 
concerns the breaking of symmetry in real space and in 
spin space. In the field-induced case the plane of the 
staggered moment is dictated by the applied magnetic 
field, and its direction within the plane is determined 
by a spontaneous breaking of the 0(2) symmetry. Ro- 
tation of the moment corresponds to the single Gold- 
stone mode. In the pressure-induced case there is no 
coupling of the spin direction to the crystal axes (in the 
absence of additional terms in the Hamiltonian which are 
not introduced here), the development of AF order re- 
quires a spontaneous breaking of 0(3) symmetry, and two 
Goldstone modes specify the rotations of the staggered 
moment. Here we have for convenience chosen the x- 
direction, but the results are identical to bond-operator 
treatments choosing any other symmetry-breaking axis. 
The second difference concerns the triplet modes partic- 
ipating in the condensate in the bond-operator descrip- 
tion: in the field-induced case the applied magnetic field 
lowers the symmetry of the system in a manner such 
that only one mode is driven to become massless at the 
QPT. It is this mode whose mixing into the singlet deter- 
mines both the primary ground-state properties of the 
ordered phase and the nature of the phase fluctuation 
(the Goldstone mode), while the small admixture of the 
highest triplet is required primarily for consistency. In 
the pressure-induced case the degeneracy of the magnons 
is lifted only by the spontaneous breaking of symmetry, 
and the condensate is composed only of the singlet and 
the triplet (tx) which constitutes the amplitude mode, 
i.e. the Goldstone modes (spin waves) explicitly do not 
form a part of the condensate. These two features are re- 
sponsible for the differences between pressure- and field- 
induced QPTs at the level of quantities such as critical 
exponents (below). Because the field acts to break the 
symmetries responsible for the properties of the pressure- 
induced QPT, any studies of quantum criticality under 
simultaneous application of finite field and pressure are 
dominated by the behavior of the field-induced QPT. 
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FIG. 14: Phase diagram for the model in the P-H plane. 

B. Finite magnetic field 

Magnetic fields and hydrostatic pressure may be re- 
garded as cooperating factors in the stabilization of an 
AF ordered phase. One may thus determine the phase 
boundary between magnetically ordered and disordered 
phases, and this is shown in the estimated P-H phase 
diagram in Fig. 1141 The phase boundary is determined 
by the condition 

g^^BH = ^J^~J[Ji+piP){J2 + 2J3)]. (60) 

The leading linear dependence of the superexchange on 
pressure in the vicinity of Pc results in the phase bound- 
ary having the form He {Pc — PY^^ ■ We note that 
the phases on both sides of the boundary show consid- 
erable evolution across the phase diagram in properties 
such as the uniform and staggered magnetic moments, 
and the splitting and dispersion of the magnon modes, 
as a function of the ratio of H to P . 

As noted in the preceding subsection, the degeneracy 
and dispersion of the magnon modes in the pressure- 
induced ordered phase are altered significantly in the 
presence of a magnetic field. A finite field lifts the de- 
generacy of the transverse spin waves because their spin 
fluctuation directions (y and z here) are no longer equiva- 
lent. The field causes the mixing of the amplitude mode, 
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FIG. 15: Magnon dispersion relations in the pressure-induced 
AF pliase under a magnetic field H = 5T < He- 
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FIG. 16: Magnetic-field dependence of spin gaps in the 
pressure-induced AF phase, illustrated for p = 1.05. The 
dashed line shows the behavior of the highest mode at the 
critical pressure. The spin gap of the middle branch, which is 
a spin wave at /i = 0, is the same for p = and p > pc- 

E]f_^^, with the phase mode i?fc_, and the mixing coeffi- 
cient changes from = 7r/4 at zero field to (/i = in the 
high-field limit, where the modes t+ and t- are again 
the appropriate eigenstates. Fig. [T^ shows the dispersion 
relations for the three nondegenerate modes at a pres- 
sure p > Pc for a finite magnetic field: the qualitative 
features of the spectra are manifestly dominated by the 
applied field, and only one Goldstone mode remains. In 
Fig. 1151 we have redefined the lower and upper magnon 
branches as E^._^ and respectively, as in the study 

of field-induced order (Sec. III). 

Fig. [THI shows the field dependence of the spin gaps of 
the three branches in the pressure-induced ordered phase 
p > Pc, with the same labeling convention as before. The 
middle branch does not mix with the upper and lower 
branches even at a finite field, so Eg^ shows a linear in- 
crease from zero which is in fact independent of p for all 
p > Pc- Because the amplitude mode has a finite ex- 
citation gap for p > Pc aX h = 0, the gap grows only 
quadratically at small fields (linearly at p = Pc)- It is 
clear that the amplitude mode is adiabatically connected 
to the upper magnon branch of the field-only problem 
(p < Pc), which for p > pc is retrieved in the high-field 
limit by a continuous alteration of the mixing coefficient 
(j) under the external magnetic field. 

C. Magnetization 

We conclude this section by discussing the uniform and 
staggered magnetization of TlCuCla under finite pres- 
sure, for which we define the magnetization exponent 
near the critical field by 

SA4 = M{H)-M{Hc)oc{H-Hcr. (61) 

Fig.llTlshows the behavior of both magnetization compo- 
nents for the three regimes p < Pc, p = Pc, and p > pc- In 
the disordered phase [Fig. llTf a^]. the critical field hc{p) is 
finite [Fig. I14| and the uniform magnetization Af^ shows 
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FIG. 17: Normalized magnetization curves for TlCuCls: (a) 
in the disordered phase (p = 1 < pc); (b) at the critical 
pressure (p — Pc); (c) in the pressure-induced AF-ordered 
phase {p — 1.05 > Pc)- Mz and Mxy denote respectively the 
uniform and staggered magnetizations, and we have taken a 
g-factor of 2.16. 





P <Pc 


P=Pc 


P> Pc 




7 = 1 
7=^ 


7 = 3 
7 = 1 


7 = 1 
7 = 2 



TABLE II: Exponent for magnetization curves. 



in the disordered phase. Increasing the pressure causes 
a decrease in Hc{p) until it reaches zero at the criti- 
cal pressure [HciPc) = 0] as in Fig. ^1 In this case 
[Fig. llZt b)] both magnetizations develop from H — 0: 
the field induces a linear mixing of singlets and triplets 
in the ground state, v (x H, and Mz depends also on the 
difference — of the weights of the t+ and t_ modes 
[Eq. (|43|) ]. whose field-induced splitting is linear whence 
7 = 3. Mxy is not sensitive to the latter splitting [Eq. 
(|44|l ] . and depends only linearly on the singlet-triplet 
admixture, yielding 7=1. Finally, the ordered phase at 
p > Pc [Fig. [rTT cl] has a finite staggered magnetization 
at zero field, and has a finite value ai H = 0, lead- 
ing to 7 = 1 for Mxy The uniform component appears 
to show a similar deviation from linearity as in the dis- 
ordered phase, but for similar reasons the leading field 
dependence is in fact linear (7 — 1) close to H = 0. 

The behavior of the two magnetization components in 
the vicinity of the critical fields is summarized by 



SMz 



J'(3J-4J')(^ ^c) 



(2J')^ 



(P < Pc) 
{P = Pc) 
{P > Pc) 



(62) 



for the uniform component, and 
1/2 



SMx 



Vf^-2J.J' 
J'(3J-4./') 

1 h 

3./-2J' 



V2J 



{p<pc) 

{P = Pc) 
h"^ {p > Pa) 



(63) 



for the staggered component, where J' denotes the quan- 
tity \^q\ defined in Eq. H54|) . The magnetization expo- 
nents are collected in Table II. 



SUMMARY 



an apparent nonlinear field dependence over a significant 
field range, as discussed in Sec. IIIB. However, close to 
the critical field the leading power-law dependence is lin- 
ear (7 = 1), reflecting the fact that the difference between 
TlCuCls and KCuCla arising from stronger interdimer 
coupling is only quantitative in nature. The staggered 
moment, M^y, has the mean-field exponent 7 = 1/2 ex- 
pected for a gapped phase (Sec. IIIB). These exponents 
may also be deduced by considering the dependence of 
the magnetizations on the coefficient [Eqs. (|43|l and 
(|44|l ]. which specifies the mixing of the triplet magnons 
into the singlet ground state (|23|l . given that cx H—Hc 



We have used a bond-operator model for TlCuCla at 
zero temperature to provide a complete description of 
quantum phase transitions (QPTs) in interacting spin- 
dimer systems driven by both magnetic field and hy- 
drostatic pressure. The magnetically ordered and disor- 
dered phases were characterized by analyzing both the 
dispersion relations of the triplet magnon modes and 
the staggered and uniform magnetizations as functions 
of field and pressure. At zero field and ambient pres- 
sure, a minimal coupled-dimer model gives a very good 
account of the magnon dispersions both in TlCuCls, de- 
spite 3d interdimer interactions which are almost suffi- 
ciently strong to close the spin gap, and in the structural 
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analog KCuCla where inter dimer coupling is weaker. At 
finite fields the bond-operator framework provides a con- 
tinuous description of the low-field, dimer-singlet phase, 
the intermediate-field regime of partial (staggered and 
uniform) magnetic order, and the fully polarized, classi- 
cal phase at high fields. Under pressure the spin gap of 
the quantum dimer phase decreases monotonically to a 
single transition into a phase of staggered magnetic or- 
der. The evolution of the order parameters through the 
QCPs separating these regimes is continuous. 

Considering first the gapped system at ambient pres- 
sure, the magnetic excitations at low fields are merely the 
Zeeman-split magnons of the zero-field case, whose dis- 
persion is unaltered by any field up to the lower critical 
field He, where the spin gap is closed. In the high-field 
regime above H^, the excitations are precisely the gapped 
spin waves expected in a fully polarized antiferromag- 
net. At He, a QPT takes place to an intermediate-field 
phase characterized by an antiferromagnetic order of the 
transverse component of the magnetization, and by one 
excitation mode which is massless at the Bragg point. 
This magnon has a ready interpretation as the phase 
fluctuation mode (Goldstone mode) of a condensate of 
mixed singlet and triplet character on each dimer rung 
for fields H^ < H < Hg. Because of its gapless nature 
the phase mode dominates the response of the system at 
intermediate magnetic fields, including most notably the 
field-dependence of the magnetization close to the QCPs. 
The bond-operator description of the ground state is a 
linear combination of singlet and triplets which changes 
continuously from pure singlet character at He to a pure 
triplet state at Hg. The description of the phase mode 
is in good quantitative agreement with recent INS mea- 
surements performed at iJ > Hc^ 

If the QPT is interpreted as a destabilization of the dis- 
ordered phase when one or more of its magnon modes be- 
come soft, a second route to an ordered phase is provided 
by the application of hydrostatic pressure to enhance the 
interdimer interactions and broaden the magnon band. 
In this situation the gapped magnons of the quantum 
dimer phase remain triply degenerate as the gap de- 
creases, until at the quantum critical point one obtains 
a phase with three spin waves. A further increase of 
the interdimer couplings leads to a conventional anti- 
ferromagnet with two spin waves - gapless (Goldstone), 
transverse, phase modes of the staggered order parame- 
ter - accompanied by one longitudinal mode which de- 
velops a gap and corresponds to amplitude fluctuations 
of the staggered moment. The appearance of pressure- 
induced order in TlCuCla has recently been confirmed 
by experiments^'' which enable one to conclude that the 
critical pressure for the QPT is very small (a result now 
confirmed experimentally-'^^), thus presenting a valuable 
additional control parameter for the experimental char- 
acterization of QPTs in the TlCuCla system. 

Both magnetic field- and pressure-induced ordering 
transitions can be interpreted as a condensation of low- 
lying magnon modes, whose linear combination with the 



singlet describes the ground state of the ordered phase. 
Indeed, the fact that the modes become soft at wave vec- 
tor Q = (0, 0, 27r) in both cases suggests a strong simi- 
larity between the QPTs driven by both magnetic field 
and pressure. However, we have shown also that an im- 
portant, if quantitative, difference exists concerning the 
degeneracy of the magnon modes and the breaking of 
rotational symmetry in the ordered phase. While the 
pressure-induced transition is characterized by two Gold- 
stone modes and a spontaneous breaking of the full 0(3) 
symmetry of spin space, the field-induced transition has 
only one Goldstone mode corresponding to a spontaneous 
breaking of 0{2) symmetry in the plane normal to the 
applied field. 

Returning to the description of the ordered mag- 
netic phase as a Bose-Einstein condensation (BEC) of 
magnons, the bond-operator formulation provides con- 
siderable additional insight into this approach. The 
Goldstone mode of the ordered phase corresponds to the 
linearly dispersive "sound" mode connected with BEC. 
The associated creation and annihilation processes arise 
from the triplet -triplet interaction Titt, which in the 
ordered phase generates terms of the form a^&fc-|-&_fc+ 
[Eq. H36|l ] as a consequence of the condensation of the 
mode t-|_. A recent numerical analysis of the scaling ex- 
ponents of the magnetization in the ordered regimc^'^ has 
confirmed that the BEC scenario provides a suitable de- 
scription of the QPT, and is universally valid for all in- 
terdimer couplings. We emphasize, however, that away 
from the transition it is necessary to modify the ground- 
state condensate, to a linear combination of the dimer 
singlet and two of the triplet components, in order to 
describe the ordered phase across the entire parameter 
region. 

The theoretical analysis of QPTs is often limited by 
the absence of suitable descriptions valid simultaneously 
for the phases on both sides of the transition. Approx- 
imate approaches applied in one phase or the other are 
fated not to converge on a single transition point, and 
thus cannot provide a continuous evolution in physical 
properties across the transition. The bond-operator for- 
mulation for the magnetic QPTs in spin-dimer systems is 
appropriate for both disordered and ordered phases, and 
thus represents an exceptional case illustrating the conti- 
nuity of physical properties across such transitions. Here 
we note, however, in the context of TlCuCls, that two 
recent, high-resolution experimental studies^^L^ have re- 
vealed a small discontinuity at the ordering transition. 
The weakly first-order magnetic phase transition is ac- 
companied by a simultaneous structural distortion, pre- 
sumably as a consequence of magnetoelastic coupling 
which is not included in the present analysis. The na- 
ture of this coupled transition awaits further investiga- 
tion. We comment briefly that the remarkable success 
of the bond-operator framework for the global phase di- 
agrams of TlCuCls and KCuCls rests to a considerable 
extent on the fact that the dimer system has significant 
coupling in all three spatial dimensions, as a result of 
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which the mean-field approximations employed here are 
highly accurate. This is reflected in the large value of the 
singlet condensate parameters in the disordered phases of 
TlCuCla and KCuCla, which deviate from unity due to 
quantum fluctuations by at most 3%. The dimensional- 
ity of the system is also responsible for the success of the 
approximation that the constraint is applied only glob- 
ally. This approach remains appropriate if it is possible 
for dimer triplets, which are effectively hard-core bosons, 
to avoid each other by perpendicular hopping processes. 
Only in one dimension (Id) is this no longer possible, and 
so the bond-operator treatment of the isolated spin lad- 
der would require a more systematic application of the 
local constraint in order to reproduce the incommensu- 
rate excitations unique to Id. 

We close by summarizing the experimental possibili- 
ties offered by the study of QPTs in magnetic systems. 
While some experiments are now under way to investi- 
gate the field-induced QPT, most notably by inelastic 
neutron scattering to measure the magnon modes both 
below and above the transition, much scope remains for 
deeper analysis by techniques including NMR. The appli- 
cation of pressure offers a very promising field for future 
studies. The experimental result of robust magnetic or- 
der at high pressure shows that it is possible to achieve a 
genuine tuning of superexchange constants through the 
QPT by the application of hydrostatic pressure. Fur- 
ther, the rather small pressures required imply the pos- 
sibility of a variety measurements, including pressure- 
dependent specific heat, magnetization, magnetic suscep- 
tibility, elastic neutron scattering for the magnetic order 
parameter and inelastic neutron scattering for magnon 
dispersion relations. All of these quantities may be mea- 
sured through the continuous QPT, and in combination 
with ongoing field-dependent studies yield the possibility 
of obtaining a complete experimental map of the QPT 
from magnetic order to quantum disorder in TlCuCls. 
We have provided a number of predictions for the prop- 
erties of the phases in this map, including most notably 
the evolution of the magnon dispersion relations and the 
magnetization exponents furnished in Table II. In this 
connection we also stress the possibility of performing ex- 
periments to detect the amplitude mode in the pressure- 
induced ordered phase, which should be visible in sus- 
ceptibility and NMR measurements, as well as by inelas- 
tic neutron scattering. In conclusion, the phenomena 



of field- and pressure-induced magnetic order provide 
excellent examples of quantum phase transitions readily 
controlled by experiment. 
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APPENDIX 

Here we present the expressions for e^.j and Af^^ used 
in the transformed Hamiltonian l|26|) at quadratic order 
in the triplet modes (|36|l . The energies are given by 

sfc- = e+g^ + e^f-v^ie+f + e^g^) (A.l) 
+v\f ~ g^ feQ + Co + - 4.fg^v^)ek, 

efc± = + e-)/.9 

-uv\l-2fg)if-g')iek + eQ), 
ta = J — ah, (q! = +,0,— ) 

where is given by Eq. (|19|l and 

Co = -eQ[2u\\\ + 2fg) - v\f - .g^)^], (A.2) 
and the coefficients of the pairing terms by 

Afc+ = -[2u\'^2fgiu* + v*) + u'v'{f~gYH, 

Afco = (u^ + 2fgv^)ek, 

Afc- = -{u^ + v^2fg)2fgek, (A.3) 

Afc± = {u' + 2fguv'){f-g')ek. 
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